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For a one dimensional half-infinite chain with diagonal disorder we calculated the 
ultimate at t — > oo value of the average excitation density at the edge site D if at 
t = the excitation was localised at the edge site (Anderson' s creterium). We 
obtained the following results: i) for the binary disordered chain we derived the 
close expression for D which is exact in the limit of low concentration of defects and 
is valid for an arbitrary energy of defects e. In this case D demonstrated the non 
analytical dependence on e. ii) The close expression for D is obtained for the case 
of an arbitrary small disorder, hi) The relative contribution of states with specified 
energy to D is calculated. All the results obtained are in complete agreement with 
computer simulation. 



I. INTRODUCTION, THE PROBLEM SETTING, AND MAIN RESULTS 

The mathematical models of contemporary physics of disordered systems can be divided 
into two classes - continuous and discrete. The continuous models are those in which the 
Shredinger equation (—A + v{r))ip = Eip with random potential v(r) is studied while the 
discrete models are dealing with the random matrix of Hamiltonian. These two types of 
models have much in common but despite the similarity they may require essentially different 
methods of analysis due to the following reason. It is well known |l( that the possibility of 
the localized states whose wave functions are essentially differes from zero within the finite 
region when the system volume runs to infinity is one of the most important properties of the 
homogeneous disordered systems. For the continuous models of infinite volume the localized 
states can be recognized by dividing the energy spectrum into discrete and continuous parts 
with square integrated states of the discrete part considered to be localized. The energy 
spectrum of discrete models being the spectrum of random matrix is always discrete and for 
this reason one should use another criterium for states characterization in this case. The 
Anderson criterium is one of the relevant onesp], 12]. 

The deep theoretical analysis is possible for the one dimensional models of disordered 
systems which we have in mind hereafter. It is commonly accepted that the spectrum of 
one dimensional Shredinger equation with random potential is completely discrete what is 
correspond to localization of all states for an arbitrary low disorder [lj. The validity of 
this statement to the discrete models is arguable because of completely discrete character 
of spectrum of such models even in the absence of disorder when all states are delocalized. 
In present paper we study the character of states in the sense of Anderson criterium for the 
case of one dimensional chain (discrete model) with diagonal disorder. The simple physical 
sense of this criterium allows one to apply it both to discrete and continuous models. 

Let's turn to the problem setting and consider the one dimensional diagonally disordered 
discrete model with the Hamiltonian whose matrix has the following entries: 



H rr > = 5 r ye r + 5 r y + i + 8 r y^i, r, r' = 1, N (1) 
Such Hamiltonian describe the Frenkel exciton in the chain consisting of N two-level 



2 



atoms in the nearest neighbour approximation. In this model the level splitting e r of r-th 
atom is considered to be the random variable distributed in accordance with the function P 
which is supposed known. 

The unit non diagonal elements of Eq.(l) are specify the scale of energy. Everywhere 
below the thermodynamic limit iV — > oo is implied. For this model we now set the following 
problem. Let us suppose that the edge atom was prepared in the excited state at t — 
and one should calculate the probability D for this atom to remain in the excited state at 
t — > oo. From the mathematical point of view it means that the initial state of the system 
is described by the wavefunction (column- vector) \£(0) with entries \l/ r (0) = 5 r ^ and one 
should find D = (\^/^(t — > oo)| 2 ) where angular brackets denote the averaging over the 
random level splittings e r . The temporary dependance of the wavefunction can be written 

as: iff(t) = exp ^Ht^^(O). Consequently the quantity of interest D can be expressed in 

terms of eigen vectors and eigen numbers Ex, A = 1, N of the matrix Eq.(l) as: 




The similar quantities were analysed in the numerical study of disordered chains for 
instance in j3[. Regarding D one can make the following qualitative conclusions. Suppose 
all the eigen functions of Hamiltonian Eq.(l) are delocalized in the sense of their amplitude 
being nearly the same within an arbitrary region of the chain. The amplitude squared of 
such functions at the edge site can be estimated as \^^\ 2 ~ 1/N with all N eigen functions 
giving nearly the same contribution \^^\ 4 ~ 1/N 2 to Eq.(2). Consequently, if the states of 
Eq.(l) are delocalized in the above sense then in the thermodynamic limit D ~ 1/N — > 0. 
Let's consider now the situation when some of the eigen functions of Eq.(l) are localized in 
the sense of their amplitude being essentially non-zero only in some restricted region of the 
chain with the size of this region being independent on N when N oo. Being defined only 
by the functions whose amplitude is essentially non-zero at the edge site the contribution of 
functions of this type to Eq.(2) will not depend on N. In this case D remain finite in the 
thermodynamic limit. Having in mind all above qualitative conclusions one can introduce 
Anderson's criterium as: If D is remain finite in the thermodynamic limit then in the set of 
eigen functions of Hamiltonian Eq. (1) there exist ones localized in the sense of Anderson's 
criterium. 

To get some information concerning the degree of localization of eigen vectors of Eq.(l) 
in some (specified) spectral range we now introduce the "participation" function defined as: 

W(U)dU=( £ |^|A (3) 

Obviously D = f W{U)dU . Reasoning analogous to that presented above shows that if 
all the states within the interval [U, U + dU] are delocalized then W(U) = 0. In the opposite 
case W(U) differs from zero. Moreover, the function (3) provide quantitative information 
concerning the average value of eigen vectors of random matrix Eq.(l) at the edge site within 
the spectral interval [U, U + dU\. The fantastic capabilities of modern personal computers 
allows one to perform direct diagonalisation of the matrix Eq.(l) for iV ~ 1000 and more 
and in such numerical experiment to "observe" quantities Eq.(2) and Eq.(3). The theoretical 
calculation of these quantities is the main goal of this paper. 
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The main results obtained in the present paper are: 

1. The perturbation theory for calculation of joint distribution function of advanced and 
retarded Green's functions of Hamiltonian Eq.(l) is developed. 

2. For the binary disordered one dimensional system, described by Hamiltonian (1) with 
the atomic splittings e r being equal to zero with probability 1 — c and equal to e with 
probability c (0 < c < 1) the following expressions for D and W(U) function are obtained: 

d= £ ij u{4 - 1,2)3/2 Ki^p +i ) +c e<|£| - i] + ° {c2) w 

W(U) = ^0(2 - |C/|)(4 - U 2 f' 2 ln(^^ + l) + (5) 

+cQ(\e\ - 1) (^)^ (U - ^±!) + 0(c 2 ) 

Nonanalyticity with respect to e is related to the occurrence (when \e\ > 1) of the edge state 
with eigen energy U = (e 2 + l)/s. 

3. For the wide class of disordered systems described by Hamiltonian (1) with random 
atomic splittings e r having the distribution functions in the form Pa{x) = p(x/A)/A (where 
p(x) > 0, J p(x)dx = 1, J p(x)xdx = 0, / p(x)x 2 dx = M 2 ) the following expressions for D and 
W(U) function are obtained: 

A 2 Jul A 2 A/f 

D = -—^ + 0(A 3 ), W(U) = 0(2 - \U\)^j^VI^TP + 0(A 3 ) (6) 
2 An 

All the results obtained in the paper were verified by the computer simulation which 
shows that when D < 0.1 the deviation of formulas Eq.(4 - 6) from the numerical results is 
less than 10%. In this case the degree of disorder may be rather large. For example formula 
Eq.(6) is still valid when e r values are uniformly distributed within the interval [—0.25, 0.25]. 



II. STATISTICS OF GREEN'S FUNCTIONS 

It is easy to show that the quantity (\^/ N (t)\ 2 ) (with the quantity of interest D = 
{ |\J/jv(oo) | 2 ) ) can be calculated as: 

(|*jv(*)| 2 ) = lim [ dU 1 dU 2 exp[t(U 1 -U 2 )t}( 1 (U 1 -tV 1 ) 1 (U 2 + tV 2 )} (7) 

Where 7(fi) - is the edge Green's function (EGF) for Hamiltonian Eq.(l): 

lv|> A I 2 

To calculate the mean product of two Green's functions entering Eq.(7) one should know 
their joint distribution function. We obtain the equation for this function by generalising 
the Dyson's method [H, 0]. Let's denote by 7(f2j),z = 1,2 the edge Green's functions Eq.(8) 
with complex energies Q± = U\ — iV\ and Q 2 = U 2 + iV 2 and add one more atom with 
splitting e to the chain. Then as it is shown in [l|, |j] EGF of the chain with added atom 
7(fi) can be expressed in terms of EGF of the initial chain as: 
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* n > = n-.- 7 (n) (9) 

To describe the EGFs of the initial chain 7(^3), « = 1,2 we introduce the distribution 
function 77 defined in such a way that the quantity 77(^1, y±, x 2 , y 2 )dx 1 dy 1 dx 2 dy 2 gives the 
probability of Re 7(^3) G [xi,Xi + dxi] and Im 7(0^ e [j/*, t/i + = 1,2. We denote by 

77 the analagous function for the chain with added atom. The relation Eq.(9) allows one to 
express fj in terms of 77 and the distribution function of atomic splittings P(e): 

rj(xi,yi,x 2 ,y2)= [ tf(Wi-Re- Wyi - Im- ) x (10) 

J V sZi — £ — xi — iy\J \ \l\ — e — x\ — iy\J 

s(x 2 - Re- V(Vs - Im- S )r](x 1 y 1 x 2 y 2 )P(e)dx 1 dy 1 dx 2 dy 2 de 

\ ll 2 - e - x 2 - %y 2 ) \ \l 2 - e - x 2 - %y 2 ) 

In the thermodynamic limit N — > 00 it must be 77 = fj. Calculating the integrals with 
(^-functions in Eq.(10) we obtain for the steady state function 77 the following equation: 

(x\ + y\f{x\ + yi) 2 r](x u 7/1, x 2 , y 2 ) = (11) 

Using the function 77 one can introduce the mean product of the advanced and retarded 
Green's functions entering Eq.(7) as a sum of four terms: 

(7(^1)7(^2)) = J i](x 1 y 1 x 2 y 2 )[x 1 x 2 - y x y 2 + %{y x x 2 + y 2 x 1 )]dx 1 dx 2 dy 1 dy 2 = (12) 

= (xix 2 ) - {y x y 2 ) + i{yix 2 ) + i{y 2 x 1 ) 

Using the expression Eq.(8) it is easy to see that these terms can be written in the 
representation of Hamiltonian Eq.(l) as: 

{ } " \ h m - e x ? + v?\[{u 2 - Ex ,y + vi] / [U} 

\ t? K^i " Ex) 2 + V?][(tf a - £y) 2 + V 2 2 ] 
f v /y g^ggM^gv) 

n^ 2 ; «\2. [(c/l _ Sa)2 + n2][([/2 _ Ea , )2 + vi] 

l**l 2 Wl 2 Wi - s A ) 



.(X ltt! > [(f/i _ + ^2 ][([/2 _ Sa , )2 + 

Thus, (|\l/Ar(/;)| 2 ) Eq.(7) can be represented as a sum of four contributions: 

(l^(t)l 2 ) = A {X1X2) + A {yiy2) + A (yiX2> + A (X1W> . (14) 
We present the calculation of A^ 2Xl ) as an example: 
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A( y2 xi) = ~r~5 li m / dU l dU 2 expt(Ui - U 2 )t (y 2 ^i) = (15) 

47T Z Vi,2—*+0 J 

lim / dU,dU 2 exptiU, - U 2 )t I £ 



' x .\ 2 \H\ 2 V2(Ui-Ex) 



AA' 

1 1 ' -,A |2| lT ,A'|2 . 



= ^E^nI^nI 2 ^^- Ey)t 



AA' 



The similar calculations shows that all four contributions Eq.(14) to (|\l/jv(t)| 2 ) are equal 
to each other: A< X1X2 ) = \ yiy2 ) = \ yiX2 ) = \x iy2 ) and therefore: 

{\V N (t)\ 2 )=4A {y2Xl) (16) 

The following remark is important for what is hearafter. Suppose that the integration 
over U\ t 2 in Eq.(15) runs over the small region C/ lj2 e [U, U + dU] only and we are interested 
in the behaviour of A^ 2Xl ) contribution at t — > oo. In this case the sum in the last string 
of Eq.(15) will contain the states with energies E\ e [£/,[/ + c#7] only. So, it is seen that 
such restriction of the integration region allows one to calculate the "participation" function 
Eq.(3). 

Thus, the problem reduced to solving of the Eq.(ll) for the joint probability distribution 
function rj. The fact that when calculating the contributions Eq.(13) the limit Vi,2 — > +0 is 
implied one can use for reducing the problem to studying the equation which is much easier 
than Eq.(ll). To do this we note that if V 1 = V 2 = then the solution of Eq.(ll) can be 
presented in the form: 



= S(y 1 )S(y 2 )p(x 1 ,x 2 ), (17) 

Vl,2=0 



r](x 1 ,y 1 ,x 2 ,y 2 ) 

where the depending on U^ 2 function p(x±,x 2 ) satisfy the following equation: 

xjx 2 2 p(x 1 , x 2 ) = J P{e)p{U 1 - e- l/x 1 , U 2 -e- l/x 2 )de (18) 

Now let us perform the calculation of the quantity (y 2 Xi) taking into account that for 
small Vi,2 the solution of Eq.(ll) goes to Eq.(17). Using the fact that the function rj is 
satisfy to Eq.(ll) one can write the following expression for the mean of interest (y 2 Xi): 

(V2Xi) = J r}{x\yxx 2 y 2 )y 2 xxdx x dx 2 dyxdy 2 = (19) 

dx x dx 2 dy x dy 2 dey 2 x x f _ _ x x y x _ x 2 y 2 \ 

(xl + yl) 2 (x 2 2 + y 2 2) 2V \ 1 £ x\ + yV Vl + x 2 + yV U2 £ xl + yV V2+ x 2 + y 2 ) n£} - 
By replacing the variables: 

2~r~2 x ^ TT 2 -> Vk k=l,2 (20) 

Vk + X k Vk + X k 

and calculating the corresponding Jacobians one can continue the equality Eq.(19) as: 

= - / dx x dx 2 dy x dy 2 de 77(1/1 - e + an,yi - Vi, E/2 - e + z 2 , V 2 + y 2 ) o^ 1 2 2 2 2 

2/i + x{ yi + xj 
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Taking advantage of the fact that when V 12 — > +0 the function r\ is close to Eq.(17) 
one can conclude that in the region where the expression under this integral is essentially 
differs from zero the following estimations are valid: y± ~ Vi,y 2 ~ —V2, with the accuracy 
of estimations increasing in the limit 14,2 - ► +0. In this limit 2/2 / (2/2 + ^1) — ¥ ~ n ^( x 2), an d 
yi — > +0. Having all this in mind one can now perform the integration over x 2 , replace the 
function 77 by its ultimate expression Eq.(17) and finally get: 

/dx 
de— p(C/i -e + x,U 2 - e)P(e), (21) 
Ob 

Here we imply the main value of the integral. The same calculations can be performed for 
(2/1^2), (2/12/2) and (x!X 2 ). 

In the end one should take into account the following important remark. It is seen from 
the Eq.(15) that if (y 2 Xi) has no singularity at U 2 — U± = uj — then lim^oo A^ 2Xl ) = 0. 
Therefore the non-zero value of D is related to the occuarence of the singularity of (y 2 xi) at 
uj = 0. Thus, for calculation of D it is sufficient to solve the Eq.(18) for small uj, extract the 
singular part and use it for the calculation of ultimate at t — > 00 behaviour of the integral 
in Eq.(15). This will be done in the next sections 



III. BINARY DISORDER 

In the case of binary disorder mentioned in the first section the distribution function of 
atomic levels splittings has the form: 

P(y) = (l-c)6(y) + c5(y-e), 0<c<l (22) 
and Eq.(18) can be written as: 

x\x\p(xix 2 ) = (1 - c)p(U 1 - 1/xi, U 2 - l/x 2 ) + cp(Ui - e - 1/xi, U 2 - e - l/x 2 ) (23) 

We now introduce the function p as a series in powers of c: 

00 



P = E <?Pn (24) 



n=0 



By substituting this series into Eq.(23) and equating the coefficients at equal powers of c 
one can obtain: 

c° : p (xix 2 )xlxl = p (Ui - l/x 1 , U 2 - l/x 2 ) (25) 

c 1 : pi(xix 2 )x\x\ = pi{Ui - 1/xi, U 2 - l/x 2 ) + (26) 

+Po{Ui — e — 1/xi, U 2 — e — l/x 2 ) - p (Ui - l/x-i, U 2 - l/x 2 ), and so on 

To calculate the value of D up to the terms of ~ c 2 it is sufficient to calculate the mean 
Eq.(21) with the same accuracy. By substituting Eq.(22) and Eq.(24) to Eq.(21) we obtain: 

(y 2Xl ) = ir J ^{(l-c)po(U 1 +x,U 2 ) + c Po (U 1 -e+x,U 2 -e) + cp 1 (U 1 +x,U 2 )} + 0(c 2 ) (27) 

The first term in braces gives the mean {y 2 X\) (up to the factor 1 — c) for the completely 
ordered chain (when D — 0) and is not of interest for us. Thus, for the singular part of 
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(y2%i) ( we denote it as "sing") which we are interested in and for the ultimate at t — > oo 
value of ^■(y 2Xl ) w e obtain the following expressions: 

sing(y 2 zi) = C7T f-^r{po{Ui -e + x,U 2 -e) + p 1 (U 1 + x, £7 2 )| + 0(c 2 ) (28) 
lim A, y Xl) = lim / dUidU 2 exp i(Ui - t^)* sing(y 2 £i) 



To obtain the functions p an d P\ entering Eq.(28) one should solve the equations Eq.(25) 
and Eq.(26). We start the analysis of these equations from the most important case when 
< 2, i.e. when energies of both Green's functions are belong to the spectrum of bare 
Hamiltonian Eq.(l) at e r = 0. 



A. The contribution of region \Ui^\ < 2 

By direct substitution one can see that in this case the exact solution of Eq.(25) can be 
written in the close form: 

p a {x 1 x 2 ) = C Ul (x 1 )C Ua (x 2 ) 1 (29) 

where Cu{x) is Lorentzian: 

Cu(x) = ^ A ~ U2 — -4 (30) 

y ' 2tt x 2 - Ux + 1 v ; 

Let us turn now to the Eq.(26) whose solution for \Ui^\ < 2 we will construct using the 
system of special functions proposed by the author in [BJ. Below we briefly review the results 
obtained in jl]. 

Define the depending on parameter U, \U\ < 2 linear operator Tiu which acts on an 
arbitrary function in accordance with the following definition: 

Huf(x) = ^f(U-l/x) (31) 

x z 

As it is shown in j5| the eigen functions (x) and eigen numbers A n (they can be numbered 
by integer n) are defined by the following relations: 



R* -x 
R-x 



■TJ + ? ,/4 _ Tp- 
Cri.r )G"U). A„: ( U _^ //L _ u2 ) ■ \<\,\ 1 (32) 



where 



JrL — , it — , rtrL 

2 2 



The map corresponding to operator TCu plays an important role in Eqs. (25), (26) and for 
this reason we will search for the solution of these equations in the form of the expansion 
into a set of functions Eq.(32). To do this we use the rules of expansion of an arbitrary 
function f(x) in a set of functions Eq.(32) obtained in 



/(*) = £ A n aP{x) (33) 
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with coefficients A n being denned by formulas: 



Using the functions Eq.(32) one can expand pi^a^) in series: 

Pi(xix 2 ) = J2 CnmO-jj 1 ^ (x^aff (x 2 ) (35) 

|n|+|ro|^0 

Substituting the series Eq.(35) into Eq.(26) and making use of the properties of functions 
Eq.(32) one can obtain: 

E C nm a^(x 1 )aff(x 2 )[l-\ n (U 1 )\ m (U 2 )]= (36) 

|n| + |m|^0 

X\X 2 

With the help of Eq.(33) and Eq.(34) one can expand the r.h.p. of this equation in a set of 
functions Eq.(32) and obtain the following expressions for the expansion coefficients of p\. 

_ J n (Ui,e)J m (U 2 ,e) 
Unm ~ 1 - A n (C/ 1 )A.(C/ 2 ) [ ' 

where functions J n {U,e) are defined as: 

The following (based on the properties of a^\x) functions) relations are take place: 

J n (U,0) = for n ^ 0, J n (U,e) = Jl n (U,e) \J n (U,e)\ 2 = \Jl n (U,e)\ 2 (39) 
Thus, the solution of Eq.(26) for p x has the form: 

Pi (xix 2 ) = E 1 x 777TT Tf^ ^k M (40) 

|n|+|m|^0 1 - M U l)*m{. U V 

As it was mentioned above only singular at uj — U 2 — Lq ~ part of this expression is 
needed for calculation of D. It easy to see that only terms with m = —n of sum Eq.(40) 
possess the required peculiarity at uj = 0. Now we make the following replacing of symbols: 
U\ — > U, U 2 — > C/ + and write down the expression for the denominators of these terms up 
to uj 2 : 

1 - X n (U)X_ n (U + uj) = ~^=^ + 0(uj 2 ) (41) 

Taking this into account we obtain the following expression for the singular part of Eq.(40): 

iVA-U 2 \J n {U,s)\ 2 (n) 
smgp^xixs) = - E crc/ (^i)^c/ (^2) (42) 
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Let us turn now to the Eq.(28). It is easy to see that the first item under integral do not 
contribute to the final result when IL^I < 2 due to the regularity of function p Eq.(29) 
at uj = 0. For this reason only the second item (it depends on function pi) remain to be 
considered: 

f dx 

sing(?/2^i) = sing nc J —pi(U + x,U + u) — (43) 



VA-U 2 ^ \J n (U,e)\ 2 rdx (n) 

lire — <7u ( u ) / — °u (U+x) 

LUJ n _£ Q IL J X 



IC 



n 



When evaluating the expression under integral we make the following replacing of variable 
x — > —1/x and make use of the fact that the functions ajj\x) are eigen for the operator 
Tiu Eq.(31). Now we can calculate the ultimate at t — > oo behaviour of the contribution 
Eq.(15): 



lim A 

t— >oo 



c f exp tut 

hm 



|fl,2|<2 



167T 2 



dxxa^jf 1 (x) 



(44) 

The ultimate behaviour at t — ► oo of the first integral do not depend on the region of 
integration over uj: 

' "° exptx 



. . expitut 
hm / duj 



-dx = ITT 



UJ J-oo X 

(n), 



(45) 



The calculation of the first moments of functions and the integrals Eq.(38) gives: 

r-2 



dxxa^ (x) 



1 n 

2 N 



Va^u 2 , 



\MU,e)\' 



e 2 + 4-U 2 



(46) 



Thus, Eq.(44) can be evaluated as: 



lim A 



t— »oo 



{V2Xl) 



r r 2 , 00 1 

-f / dU(A - u 2 f' 2 - 

l07T 7-2 ~t n 



|c/ 12 |< 2 167r 7-2 ' n \e 2 + 4 - U 2 

And taking into account that Yl™=i q n /n = — ln(l — q), we finally obtain 



+ 0(c 2 



(47) 



lim A 



t— »oo 



(V2X\) 



|L/ 1]2 |<2 167T 



fj U ^- U * ) V^( J f_ + ^ +0{ 



(48) 



B. The contribution of region \U\2\ > 2 



Let us return to the Eq.(28) and consider the contribution of the region |C/i )2 | > 2 to the 
integrals in Eq.(28). It is shown in the Appendix that the second item in braces in Eq.(28) 
(depending on pi) do not contribute to the final result. Thus, the contribution of the region 
\U\p\ > 2 is defined by the first (depending on p ) item of Eq.(28) only. As the expression 
for EGF of the ordered chain is known [6[ the solution of Eq.(25) for the unperturbed joint 
probability distribution function p (for \Ui^\ > 2) can be easily guessed: 



Po{xix 2 ) 



fA - sign^W^f - 4 



x 2 - 



U 2 -sign(U 2 )JUl-A 



(49) 



10 



We calculate now the contribution of the first (depending on p ) item in braces (28) (we 
denote it by Ii): 



h 



ic f dxdU\dU 2 



ic f dxdU\dU 2 



x 



exp[i(Ui - U 2 )t\ p (Ui - e + x, U 2 - e) 



An 



x 



expMUi - U 2 )t\ 5 



x — e + 



f/ 1 + sign(f/ 1 ) v /f/i 2 -4 



x 



(50) 



xS 



U 2 + sign(U 2 )JUl-A 



— e 



The integrations over U± >2 run over the regions \Ui j2 \ > 2. Integrating Eq.(50) over x we get: 



dUidU 2 expi{U 1 -U 2 )t 



ic r 



U 2 + sign(U 2 )Ju 2 2 -4 



(51) 



It is clear that the result do not depend on the sign of e and below we perform calculations 
for e > 0. In this case 5-function under the integral Eq.(51) gives zero for U 2 < —2 
and, moreover, if < e < 1 this 5-function is equal to zero identically. Consequently the 
contribution under consideration is differ from zero only when \e\ > 1. The function under 
integral in Eq.(51) has a ^-peculiarity with respect to U 2 and a pole-peculiarity with respect 
to Ui To calculate these integrals (at t — > oo) we take into account the following property 
of the relevant function $(£/) defined as: 



HU) = e 



This function is equal to zero at: 



and can be expanded at U f* Uq as: 



U = U t 



€ 2 + l 







(52) 



(53) 



>o — 



$([/) k(U-U { 

Using this formula one can perform the integration over U 2 in Eq.(51): 

-u 2 + Ju. 



J dU 2 exp[-tU 2 t}5 



q -4 



exp[— iUot] 



(54) 



(55) 



Let us now turn to the integration over U\ in Eq.(51). Due to the fact that one should 
calculate this integral at t — > oo the contribution of pole-peculiarity is the only one of 
importance. Using the expansion Eq.(54) again one can obtain the following expression for 
this integral: 



L 



dU\ exp^it] 



^> 2 e - [U x + JUl - 4]/2 



1 r+oo 
J —oo 



dUi exphUit] . TT . e 2 — 
= — 27rexp \iU t\ — 



1 



C/i-C/o 



(56) 
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The approximate equality here becomes exact at t — > oo. Thus, 

c fe 1 - 1\ 2 

^ = ^\~^—) (57) 

The additional study shows that this contribution is related to the occurrence (when \e\ > 1) 
of the edge state with the energy Uq defined by Eq.(53). To obtain the final formula for D one 
should sum the contributions Eq.(57) and Eq.(48) and in accordance with Eq.(16) multiply 
the result by four. Thus, we obtain the Eq. (4) for D for the case of binary disordered chain. 
In accordance with the logic of the above calculation and taking into account the remark 
after Eqs.(15,16) one can see that the "participation" function W(U) Eq.(3) is define by 
Eq.(5). 



IV. AN ARBITRARY SMALL DIAGONAL DISORDER 

The analysis of the binary disordered chain presented above can be regarded as a consis- 
tent perturbation theory for the statistics of advanced and retarded Green's functions with 
concentration c of atoms with the level splitting e playing the role of small parameter. The 
perturbation theory for the case of the chain with an arbitrary small diagonal disorder can 
be constructed in a similar way. The relevant small parameter can be defined as follows. 

Let the function p(e) possess the following properties: p{e) > 0, / p(e)de = 1. With the 
help of this function we construct now the following family of the atomic levels splittings 
distribution functions Pa(s)' 

= (58) 



A"VA 

If the moments of p(e) function are: 



M n = J p(e)e n de (59) 
then the moments of functions Eq.(58) can be expressed as: 

P A (e)e n de = A n M n (60) 



/ 



It is clear that one can consider A as a degree of disorder - the system becomes ordered at 
A — > 0. Without loss of generality one can say that M 1 = 0. Thus, one should construct 
the perturbation theory for the Eq.(18) when P(e) = Pa(^) with A playing the role of small 
parameter. Therefore we now construct the expansion of solution of Eq.(18) in powers of A 
and we start from the case of \Ui t2 \ < 2 which is of particular importance. For this reason 
we write the expansion of the function p(x 1 x 2 ) in the vicinity of x w ,x 20 : 



oo y Qn+m 

p(x 1 x 2 )= £ p nm (x x - x 10 ) n (x 2 - x 20 ) m , Pnm = — ———p(Q 1 Q 2 
=0 n\m\ oo 1 ou 2 



(61) 



-XiQ 



By applying this expansion to the r.h.p of Eq.(18) with Xio = Ui — l/xi,i = 1,2 one can 
obtain the expansion of the function p(U\ — e — 1/xi, U 2 — e — l/x 2 ) in powers of e and 
express the r.h.p. of Eq.(18) in terms of moments Eq.(59) and powers of A: 

dePA{e)p[Ui - e - 1/x-y, U 2 - e - l/x 2 ) = >^ - 



/ 



n,m=0 



n\m\ de\de\ 



t/i-l/xi 

(62) 



12 



It follows from this that the l.h.p. of Eq.(18) (i.e. the function p(xxx 2 ) itself) can be also 
expanded in powers of A: 

oo 

p(x 1 x 2 ) = J2A k Q k (x 1 x i ), (63) 

fc=0 

Eq.(18) allows one to express the functions Q n in terms of Q m) m < n. To do this we 
substitute the expansions Eq.(62) and Eq.(63) into Eq.(18). We get 



oo 



x\x\Y, A k Q k {xxX 2 ) = ■ 

k=0 n,m,k=0 



-l) n+m A n+m+k M n+m 8 n+m 



nimi 



Qk(9x9 2 ] 



=Ui—l/xi 



(64) 



Equating the coefficients at A we obtain: 

x\x\Qq{xxx 2 ) = Qo(Ux - 1/xx, U 2 - l/x 2 ) (65) 
And consequently (see Eq.(29)) we obtain the following expression for Qq{xxx 2 ) 

Qo{xix 2 ) = po{xix 2 ) = C Vl {xx)C U2 {x 2 ) (66) 
It is easy to see that if M\ = then Qi(xix 2 ) = 0. Equating the coefficients at A 2 we get: 



x\x\Q 2 {x x x 2 ) = Q 2 (U 1 - 1/xi, U 2 - l/x 2 ) + M 2 



d 2 lc?_ lc?_ 
89x392 + 2 86 j + 2 M\ 



Qo{9x6 2 ) 



9i=Ui-l/xi 

Let us now calculate up to the terms of ~ A 3 the value of (y 2 xx) defined by Eq.(21) with 
P(s) = Pa(e)- We use the expansions Eq.(62) and Eq.(63) for the relevant integral: 



deP^{e)p{Ux — e + x, U 2 — e) 



deP A (e) Yl 

n,m,k,=0 



(-e) n+m A k 8 n+m Q k {6x6 2 



E ■ 

n,m,k=0 



■l) n+m A n+m+k M n+m 8 n+m Q k ( 



h02 



n\m\ 



A 5 



Q 2 (Ux+x,U 2 )+M 2 



89189^ 
18 2 \8 2 



89x89 2 2 89\ 289\ 



nimi 

e 2 =u 2 



89\89\ 



e 2 =u 2 



(68) 



8 2 



Qo(Ux+x,U 2 ) + 



+ 0(A 3 





e 2 =u 2 - 


)Q {e 1 e 2 ) 


0i=Ui+x. 



The first term of zero order (with respect to A) is correspond to the chain with no 
disorder and therefore do not contribute to the value of D we are interested in. To calculate 
the contribution of the first term in square brackets to D (we call it A-term) one should 
obtain the function ^2(^1^2)- This function can be found from Eq.(67). Solution of this 
equation for \Ux j2 \ < 2 can be performed in the way similar to that for the Eq.(26) and 
reduced to the following redefinition of quantities J n (U): 



Jn{U) 



dx 8C U {9) 



x 2 G n (x) 89 



e=u-i/x 



UU) = r n (u) 



(69) 



By making the following replacing of variable in these integrals: 9 = U — 1/x and making 
use of the following property of the function G{9) Eq.(32): G(1/(U - 9)) = G{6)/X 1 Q we 
get: 

d£u(x) dx 



Jn{U) — A r 



dx G n (x) 



(70) 
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Calculation of these integrals shows that J±i(U) are the only nonzero of them: 

i\i , , I 



J±i(U) = =F 



\J±i(U)\' 



A-U 2 



(71) 



Having this in mind it easy to see that A-term from Eq.(68) (when < 2) can be written 
as: 

A 2 M 2 r 2 r-. =^.„ A 2 M 9 

A-term 



A< OT xi)(* 



oo, 



J VA^LPdU = 



(72) 



16tt y_ 2 8 

In the case of > 2 the analysis of Eq.(67) similar to that described in the Appendix 
shows that A-term is equal to zero. Consider now the contribution of the second term in 
square brackets of Eq.(68) to D (we call it S-term). When calculating this contribution the 
region Lq^ < 2 do not play any role for the ultimate (at t — > oo) behaviour of A^^} because 
of the absence of peculiarities at U\ = U 2 of the functions under integrals. For \U\ t2 \ > 2 
5-term from Eq.(68) can be written as: 



= M 2 A\[ d ^5( X + 
J x ax V 



£>-term 



U 1 + sign(U 1 )Juf-4 



x 



(73) 



d ,/ U 2 + sign(C/ 2 ) JU. 
X Ty\ V + 2 



exp i{Ui — U 2 )t 



t—xx 



+ ... 



y=o 



(Here we present only the terms with crossing derivatives. The remaining terms can be 
analyzed in the same way.) Due to the fact that the argument of the second 5-function 
never becomes zero (for |C/ 2 | > 2) we come to the conclusion that the 5-term is equal to 
zero and after multiplying the result Eq.(72) by factor of four we obtain formulas Eq.(6). 



V. NUMERICAL EXPERIMENT 

The formulas Eq.(4 - 6) can be verified by calculation of quantities D and function 
W{U) by formulas Eq.(2) and Eq.(3) where the eigen vectors \I/ A and the eigen energies E\ 
are obtained by direct computer diagonalisation of Hamiltonian Eq.(l). Below we present 
the results of such verification. The noisy curves were obtained numerically and smooth 
ones by formulas Eq.(4 - 6). Fig. la shows the dependences of quantity D/c on the energy 
of defect e obtained numerically for different concentrations of defects c and the relevant 
theoretical curve Eq.(4). It is seen that when D/c becomes independent on c the Eq.(4) is in 
complete agreement with the numerical results and the curving at e — 1 related to the non 
analytical part of Eq.(4) is well pronounced. Fig. lb shows the dependence of D on the degree 
of disorder A for the case of uniform disorder when the atomic level splitting distribution 
function has the form Eq.(58) withp(x) = 0(0.5— \x\). It is seen from this figure that Eq.(6) 
is in good agreement with numerical results even for rather strong disorder. Fig. 2 shows 
the "participation" functions W(U)dU calculated numerically by Eq.(3) and by analytical 
formulas Eq.(5) (binary disorder) and Eq.(6) (uniform disorder). Fig.2(a,b) relate to the 
case of binary disordered system with concentration of defects c = 0.03 and energy of defects 
e = 0.8 (Fig. 2a) and e = 1.2 (Fig. 2b). It is seen that for e > 1 the "participation" function 
" W(U) demonstrate the sharp maximum at U — U (Eq.(53)). Fig.2c shows the case of 
rather strong (A = 0.5) uniform disorder. It is seen that in this case a good agreement 
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between the numerical experiment and the theory is also take place but for the description 
of noticeable dip in the center of "experimental" curve one should take into account the 
corrections of higher order than A 2 . 
No fitting was performed. 



VI. APPENDIX 

To clear up the role of the second item (depending on p\ ) in Eq.(28) for \U\ t2 \ > 2 one 
should obtain the solution of Eq.(26) for this spectral region. For the sake of certainty let 
us consider the case of > 2 and introduce the following quantities: 

2 U - 7[/ - £ 

where is the EGF of the chain with no disorder. For the case under consideration ( 
U 1>2 > 2) the solution of Eq.(25) gives the following expression for the function p : 

po(x 1 x 2 ) = 5(x 1 - 7c/i)5(x 2 ~ 7t/a) ( 75 ) 
Then Eq.(26) can be rewritten as: 

Pi(xix 2 ) - ^—^ — 1 ^ X 2 1, 2 U ' 2 — ^— ^ = 5^ Xl - 1 ' u js(x2-j' U2 )-S(x 1 -j Ul )8(x 2 -^u 2 ) (76) 
The direct substitution shows that the solution of Eq.(76) has the form: 

p x [x x x 2 ) = lim V <^ 5{x 1 - 9 i (U 1 ))8{x 2 - 8i(U 2 )) - 5(xi - ^Ui)S(x 2 - 7^) > (77) 

i=l 

with quantities 6i(U) defined by the following recurrent relations: 

*i(C0=7k. n+l {U) = u _ 1 dn{uy n=l,2,..,M-l (78) 

Using Eq.(76) one can write the following expression for the function pi(U\ + x, U 2 ) entering 
Eq.(28): 

2 r 



p l (-l/x,y)-8(l/x+j' Ul )5(y-j' U2 )+S(l/x+j Ul )8{y-j U2 



(79) 



Pi(Ui+x,U 2 ) = lim - 

It is clear that the last two terms with 5-functions have zero limit at y — > 00. By substituting 
the function p\ (Eq.(77)) to this expression one can see that the limit of the first term is 
also equal to zero. Thus, the contribution of the second item in Eq.(28) is equal to zero for 
|[/i, 2 |>2. 
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FIG. 1: (a): The case of binary disorder. Noisy curves - the dependences of D/c on the energy 
of defects e obtained numerically for different concentration of defects c = 0.01,0.03,0.1, smooth 
curve - the relevant theoretical curve, (b): The case of uniform disorder. The dependence of the 
ultimate at t —* oo density of excitation D on the edge site on the degree of disorder A. 



Captures 

Fig.l (a): The case of binary disorder. Noisy curves - the dependences of D/c on the en- 
ergy of defects e obtained numerically for different concentration of defects c = 0.01, 0.03, 0.1, 
smooth curve - the relevant theoretical curve. 

(b): The case of uniform disorder. The dependence of the ultimate at t — > oo density of 
excitation D on the edge site on the degree of disorder A. 

Fig. 2 The "participation" function - the comparison of theory (smooth curves) and com- 
puter simulation (noisy curves). For the case of binary disordered chain the occurrence of 
the peculiarity related to the edge state is seen: (a) - e = 0.8 < 1, no strong peculiarity, 
(b) - e = 1.2 > 1, sharp peak appear, (c) - the "participation" function for the chain with 
uniform disorder at A = 0.5. dU = 1/50 for all cases. 
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FIG. 2: The "participation" function - the comparison of theory (smooth curves) and computer 
simulation (noisy curves). For the case of binary disordered chain the occurrence of the peculiarity 
related to the edge state is seen: (a) - e = 0.8 < 1, no strong peculiarity, (b) - e = 1.2 > 1, sharp 
peak appear, (c) - the "participation" function for the chain with uniform disorder at A = 0.5. 
dU = 1/50 for all cases. 



